Exact expression for the number of states in lattice models 
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We derive a closed-form exact combinatorial expression for the number of states in canonical 
systems with discrete energy levels. The expression results from the exact low-temperature power 
series expansion of the partition function. The approach provides interesting insights into basis of 
statistical mechanics. In particular, it is shown that in some cases the logarithm of the partition 
function may be considered the generating function for the number of internal states of energy 
clusters, which characterize system's microscopic configurations. Apart from elementary examples 
including the Poisson, geometric and negative binomial probability distributions for the energy, the 
framework is also validated against the one-dimensional Ising model. 
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It has long been recognized that the number of states, 
g{E), with a given value of the energy, E, is a key 
quantity in equilibrium statistical physics. In princi- 
ple, all energy-related thermodynamic properties of a 
classical canonical system can be calculated once g{E) 
is known. In particular, in the case of a system with 
discrete energy levels, the canonical partition function 
is Z{[i) = g{E)e~^^, where the summation is 
over the allowed energy levels, and quantities such as 
the Helmholtz free energy, F{f3) = —\nZ{/3)//3, the 
ensemble-averaged energy, {E) = —d\n Z{(3)/dl3, and 
the specific heat (which is related to energy fluctuations, 
(E^) - (E)^ = -d{E)/dl3) can be simply obtained. 

Due to the central role of the energy distribution func- 
tion, a variety of theoretical and computational studies 
have addressed the problem of how to obtain g{E), see 
e.g. [iHl]. On the other hand, relatively little was done 
to understand, how the state-space representation of a 
considered system relates to its energy distribution, and 
how to recover such an information from, for example, 
the partition function. 

In this paper we derive a new and general combinato- 
rial formula for the density of states. The formula results 
from the exact low-temperature power-series expansion 
of the canonical partition function. In our approach g{E) 
is expressed by the so-called Bell polynomials, which have 
a very convenient combinatorial interpretation: they de- 
scribe, how equal energy portions can cluster together 
within the studied system. This, in turn, gives insight 
into system's microscopic details. 

In the low-temperature expansion of the partition func- 
tion, Z{(3) = e~^^('^\ microstates, fl, are counted in 
order of their importance as the temperature is increased 
from zero 0, @| • If ground state of a system is known 
and if the successive excitations from this state can be 
classified in a simple way, one should, at least theoreti- 
cally, be able to construct such a series, i.e. 



where 



and energy is considered to be discrete 



E = sN, 



(2) 



(3) 



Z{x) = J2 9(.eN)x''^ 



(1) 



N=0 



with e representing a portion of energy and N — 
0, 1, 2, . . . (with the ground-state energy equal to zero). 

In Eq. ([T]), g{sN) represents the number of microscopic 
configurations in which energy of the system is sN. 

From Eq. ([T]) it is apparent that the exact distribu- 
tion function for the energy, g{eN), can be determined 
from the coefficients of the partition function. As a rule, 
however, when the order of expansion is increased com- 
plexity of contributing terms also increases rapidly. In 
this context, a rule of thumb states that the work in- 
volved in direct calculating the next term, (/(eiV), is the 
same as that needed to calculate all the preceding terms, 
g{0),g{e), . . . ,g{eN — e). For this reason, generation of 
lengthy low-temperature power-series for the partition 
function is though to be a highly specialized art. That is 
why, even in the case of the two-dimensional Ising model 
on a square lattice, for which the method was primary 
advanced (see ch. 6 in the exact distribution of en- 
ergies was just obtained in 1996 [2, 9], nearly a half of a 
century after Kaufman's exact solution for the partition 
function fioj. 

The aim of this paper is to show that the effort needed 
to calculate the coefficients of the low-temperature 
power-series expansion for the canonical partition func- 
tion is greatly reduced by using the combinatorial scheme 
described below. The approach is very general and ap- 
plies to any lattice model. We also argue that the frame- 
work provides interesting insights into basis of statistical 
mechanics. In particular, it is shown that in some cases 
the logarithm of the partition function, \nZ{x), may be 
considered the generating function for the number of in- 
ternal states of energy clusters, which characterize sys- 
tem's microscopic configurations. Apart from elemen- 
tary examples including the Poisson, geometric, and neg- 
ative binomial probability distributions for the energy. 
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the framework is validated against the one-dimensional 
Ising model, which is not at all a trivial example since its 
ordering temperature is zero, and the series, Eq. ([T]), is 
known to diverge at a: = as expected for a model with 
a zero temperature phase transition Q . 

Derivation of the main result. To derive the main re- 
sult of this paper, which is the general combinatorial for- 
mula for the number of states, we begin with the well- 
known relation between the canonical partition function 
and the Helmholtz free energy, i.e. 



Zix) 



(4) 



where x is the expansion variable of the low-temperature 
series, Eq. ([2]), and 



A{x) ^ lnZ{x) ^ 



(5) 



with a„ = d"^ A{x) / dx'^\x=o- Then, using the generat- 
ing function for Bell polynomials, Yj^{{an}), which are 
the polynomials with a sequence of parameters {a„} — 
oi, a2, . . . ajv, given by the formal power series expansion 
[Uii, Eq. Q can be written as 



Z{x) = exp 



CO ^ 

E ^>V({a„})x^. (6) 



JV=0 



In the last formula, the complete Bell polynomials, 
YN{{o,n}), are defined as follows 



Yo{{an}) - 1, 



and for all TV > 1 



N 



(7) 



(8) 



fc=i 



where Bpf^k({a.n}) represent the so-called partial (or in- 
complete) Bell polynomials, which can be calculated from 
the expression below 



N-k+l 



(9) 



where the summation takes place over all integers c„ > 0, 
such that 



N-k+l 



N-k+1 



Cn — k and E] — N. 



(10) 



Finally comparing Eqs. ([T]) and (|6]), one gets the exact 
expression for the number of states: 



g{eN) = —YNi{an}). 



(11) 



The basic difficulty with Eq. (jll[) may arise from an un- 
acquaintance with Bell polynomials as given by Eqs. (|S]) 



and ([9]). For this reason, we follow by explaining their 
meaning. 

Suppose that N distinguishable particles are parti- 
tioned into k non-empty and disjoint clusters of > 
elements each, where X]i=i — There are exactly 



)-^'■fl^-^^I[\-^ (12) 

' i—l n—l ^ 



of such partitions, where c„ > stands for the number of 
clusters of size n, with the largest cluster size being equal 
to iV — fc -I- 1, and where Eqs. pO)) are satisfied. Suppose 
further that in such a composition clusters of the same 
size are indistinguishable from one another, and each of 
Cn clusters of size n can be in any one of a„ > internal 
states. Then the number of partitions becomes, cf. with 
Eq. (USD, 



AT-fc + l 



TV! 



n A©' 



(13) 



Summing Eq. (jl3p over all integers c„ > specified 
by Eqs. (|10p . one gets the partial Bell polynomial, 
BN,k{{0'n}), which is defined by Eq. ([SJ. If for all 7i > 1 
the coefficients a,i > 0, the polynomial describes the 
number of partitions of a set of size N with exactly k 
subsets, where each coefficient, a„, represents the num- 
ber of internal states of a cluster of size n. Finally, sum- 
ming the partial polynomials over k one gets the complete 
Bell polynomial, ({««}), the combinatorial meaning of 
which is obvious. 

According to the explanations above the main result 
of this paper, Eq. (fTT|) . shows how equal energy portions, 
e, are distributed, and how they cluster together within 
the system whose structural details are hidden in the 
coefficients {a„}. If the coefficients are non-negative, i.e. 
Vn>ia„ > 0, they can be interpreted as thermodynamic 
probabilities of energy clusters of a given size. Otherwise, 
their meaning is not clear. 

In what follows, the Poisson and the negative binomial 
distributions are discussed as elementary mathematical 
examples of application of the combinatorial approach. 
The advantages of the new framework with respect to 
traditional methods of statistical mechanics are enhanced 
with reference to the one-dimensional Ising model. 

Poisson distribution. As the first example, the Pois- 
son distribution for the energy of a canonical system is 
considered, 



PiE) 



-{E) 



E\ 



(14) 



where, see Eq. ([3]), £ = 1 and = iV = 0, 1, 2, . . . . In 
order to bring Eq. (jl4p into the canonical form. 



P{E) = g{E) 



-PE 



(15) 
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one has to assume that {E) = e ^ . Then, the Poisson 
distribution can be written as 

1 e"^^ 

P^E) = (16) 
where (compare Eqs. ([T5|) and 

5(£;) = -l and hiZ(/3) = e-^. (17) 

To verify the combinatorial approach introduced in 
this paper, let us first notice that in the case of the 
Poisson distribution x — e~^, Eq. ([3]). Therefore, the 
logarithm of the partition function is simply, cf. Eqs. ([5]) 
and dm), 

A{x) = X. (18) 

The identity function for A{x) provides a very significant 
sequence of the coefficients {a„} for the energy distribu- 
tion function, Eq. ([TT|) . i.e. 

ao = and {a„} = 1, 0, 0, . . . , 0. (19) 

The sequence states that in the case of a canonical system 
with the Poisson distribution of energy, energy portions, 
e = 1, are independent from each other. In some sense 
it means that such a system is structureless (it may, for 
example, consists of noninteracting parts). 

Finally, inserting the obtained coefficients into Eq. pT|) 
one can easily get the energy distribution function, 

Eq. ini), 

5(A^) = ^>"jv(l,0,0,...,0) = ^, (20) 

where Y/v(l, 0, 0, . . . , 0) = 1, because there is only one 
composition of a set of size N in which only clusters of 
size n — 1 may exist (ai — 1 and a„ = for all n > 1). 

Negative binomial distribution. As the second example 
the negative binomial distribution for the energy of a 
canonical system is studied, 

P(i?)=(^+^-^)(l-x)^x^ (21) 

where x G (0, 1) and r > 0, while e = 1 and E = N = 
0,1,2,..., see Eq. (jS)). It is worth noting that among the 
discrete distributions, the negative binomial distribution, 
NB{r,x), Eq. (1^ . is considered the discrete analogue 
of the Gamma distribution, r(r, /3), whose probability 
density function is PiE) = fi'' E'-^e'^^ /r{r), where r 
is the so-called shape parameter. The analogy is even 
more apparent if one realizes that, in the case of r = 1 
the negative binomial distribution becomes the geometric 
distribution, NB{l,x), which is the discrete analogue of 
the exponential distribution, r(l,/3). 

Bringing the negative binomial distribution, Eq. (Hi}, 
into the canonical form, Eq. dTSl) . for x = e^^ one gets: 

,iE)^ (22) 



and 

y4(x) = -rln(l-x)=rV(n-l)!^. (23) 

n! 

n— 1 

The series expansion of the logarithm of the partition 
function, A{x), gives the following sequence of the coef- 
ficients {a„}: 

ao = and V„>i = r{n — 1)!. (24) 

The sequence provides a meaningful microscopic infor- 
mation about the considered canonical system. First, it 
appears that the parameter r describes the degeneracy 
of internal states of the energy-clusters, and the lack of 
degeneracy, r = 1, results in the geometric/exponential 
distribution. Second, since (n — 1)! is the number of per- 
mutations of n objects where the first object is fixed, 
the only reasonable explanation behind a„ oc (n — 1)! is 
the linear ordering of energy portions within the clusters. 
At first glance the idea of ordering of indistinguishable 
energy portions may seem unacceptable, but it gets intel- 
ligibility if one realizes that in some sense energy portions 
become distinguishable once they, for example, become 
interaction energies between given pairs of particles. 

To validate the main result of this paper, Eq. (|lip . 
one has to insert the coefficients given by Eqs. (p4)) into 
the mentioned expression. After some algebra one gets 
Eq. ([221): 

= 47>V(Mn-l)!}) (25) 



= ^E^'^^'^(0''l''2!,...) (26) 




where basic combinatorial identities have been used, in- 
cluding [HI, [l2|: i. the generating function for the un- 
signed Stirling numbers of the first kind, [^] , and ii. 
properties of Bell polynomials, i.e. B]\[^k {{fu'^an}) — 

t'^u^BnAWu}) and B^Aii'^ " 1)!}) = [fe] El- 

One- dimensional Ising model. The Hamiltonian of the 
closed chain of V Ising spins, {si}, with nearest-neighbor 
interactions in the absence of external magnetic field can 
be written as 

V 

H{{s,}) = -JY,s^s^+l, (28) 

1=1 

where J > 0, Si = ±1, and the periodic boundary con- 
dition is imposed by assuming that sy+i = si. In the 
thermodynamic limit, V ^ 1, the partition function of 
the model is given by 

Z(/3) = (e^-^ + e-^^)"^. (29) 
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To use Eq. pT|) to calculate the number of states, 
g{E), characterizing the model, one must first rescale 
its energy, E, because in the present form of the Hamil- 
tonian, Eq. (1^ . the ground state energy is negative, 
H{{-1}) = H{{+1}) = -JV, whereas it should be at 
least equal to zero, see Eqs. ([T]), ^ and ([2]). The prob- 
lem can be solved by adding to Eq. (1^ a constant, JV, 
so that for all the spin configurations, {si}, the modified 
Hamniltonian, 



do not satisfy the condition of non-negativity, which (if 
satisfied) allows to conclude on state space representa- 
tion of the considered model. Nevertheless, the knowl- 
edge of {a„} (regardless of their sign) enables a direct 
validation of the combinatorial approach. By inserting 
the coefficients into the expression for the energy distri- 
bution function, Eq. ([TT]) . after some algebra one gets, 
for N = E*\ 



H*i{s,}) = -JJ2s^s,+l + JV, 



i=l 



results in non-negative energy 

E* = E + JV >0. 



(30) 



(31) 



With this change, the partition function of the model 
becomes 



Z*(/3)=e-^-'^Z(/3) = (l + e-2/3J)V^ 



(32) 



but statistical properties of the chain of spins, {s;}, 
remain the same, as compared with the original Ising 
model. In particular, the probability distributions for 
spin configurations, P{{si}) = exp[—(3H{{si})]/Z{/3) 
and P*{{s,}) = exp[-PH*{{s^})]/Z*{|3), coincide with 
one another. The same applies to probability distri- 
butions for the energy, P{E) = g{E)e^l^^ /Z{l3) and 
P*{E*) = g*{E*)e-P^' /Z*{fi). 

In view of the above considerations one can see that 



g{E)=g*{E*)=g*{E+JV). 



(33) 



Thus, to determine the number of states for the one- 
dimensional Ising model, g{E), one can calculate g*{E*), 
and then change the variables E* and E according to 
Eq. dSI]). 

In the following, for simplicity one assumes in Eq. p2p 
that J = 1/2. Then, the logarithm of the partition func- 
tion, lnZ*(/3) can be written as, Eq. ([SJ, 



(-l)""i(n-l)!^, (34) 



where e = 1 and x = e ^ . From the series expansion, it 
is obvious that the coefficients {a„}, i.e. 

ao = and V„>ia„ = y(-l)"-i(n - 1)!, (35) 



Nl 
1 



Y^i{V{-ir-\n-l)l}) 



N 



'1) 



(36) 
(37) 



k—l 

where (— 1)^ '^[^] is the signed Stirling number of the 
first kind. From the last expression one immediately gets 



9iE) = 



V 

E + V/2 



(38) 



where E is the energy of the original Ising model de- 
scribed by Eqs. (^5]) and 



The obtained result, Eqs. (|37|) and (|38)) . is exactly the 
expected one: the number of microscopic configurations 
of the closed chain of V Ising spins is equal to the number 
of ways to choose positions for E* pairs of equal spins, 
{+1, +1} and { — 1, —1}, from the available V positions. 

Summary. Our approach differs crucially from previ- 
ous works on the number of states. We use the combina- 
torial mathematics of exponential generating functions to 
obtain the exact low-temperature power-series expansion 
of the canonical partition function. The expansion results 
in the exact expression for the energy distribution in ar- 
bitrary lattice models. Apart from elementary examples 
including the Poisson, geometric, and negative binomial 
probability distributions for the energy, the framework is 
validated against the one-dimensional Ising model. The 
work on the application of the formalism to the Ising 
model on a square lattice is in progress and will be re- 
ported in a future article to follow. 
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